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Abstract 

We analyze some systems of partial differential equations arising in the theory of mean field type 
control with congestion effects. We look for weak solutions. Our main result is the existence and 
uniqueness of suitably defined weak solutions, which are characterized as the optima of two optimal 
control problems in duality. 


1 Introduction 

In the recent years, an important research activity has been devoted to the study of stochastic differential 
games with a large number of players. In their pioneering articles [HIIIITI], J-M. Lasry and P-L. Lions 
have introduced the notion of mean field games, which describe the asymptotic behavior of stochastic 
differential games (Nash equilibria) as the number N of players tends to infinity. In these models, it is 
assumed that the agents are all identical and that an individual agent can hardly influence the outcome 
of the game. Moreover, each individual strategy is influenced by some averages of functions of the states 
of the other agents. In the limit when N +oo, a given agent feels the presence of the other agents 
through the statistical distribution of the states of the other players. Since perturbations of a single 
agent’s strategy does not influence the statistical distribution of the states, the latter acts as a parameter 
in the control problem to be solved by each agent. 

Another kind of asymptotic regime is obtained by assuming that all the agents use the same distributed 
feedback strategy and by passing to the limit as N ^ oo before optimizing the common feedback. 
Given a common feedback strategy, the asymptotics are given by the McKean-Vlasov theory, [MHO] 

: the dynamics of a given agent is found by solving a stochastic differential equation with coefficients 
depending on a mean field, namely the statistical distribution of the states, which may also affect the 
objective function. Since the feedback strategy is common to all agents, perturbations of the latter 
affect the mean field. Then, having each player optimize its objective function amounts to solving a 
control problem driven by the McKean-Vlasov dynamics. The latter is named control of McKean-Vlasov 
dynamics by R. Carmona and F. Delarue laiz] and mean field type control by A. Bensoussan et al, Hi]. 
When the dynamics of the players are independent stochastic processes, both mean field games and control 
of McKean-Vlasov dynamics naturally lead to a coupled system of partial differential equations, a forward 
Fokker-Planck equation and a backward Hamilton-Jacobi-Belhnan equation. For mean field games, the 
coupled system of partial differential equations has been studied by Lasry and Lions in [Hill [Hj. 
Besides, many important aspects of the mathematical theory developed by J-M. Lasry and P-L. Lions 
on MFG are not published in journals or books, but can be found in the videos of the lectures of P-L. 
Lions at College de France: see the web site of College de France, m- One can also see m for a brief 
survey, and we mention a very nice article on weak solutions of Fokker-Planck equations and of MFG 
systems of partial differential equations. 

The analysis of the system of partial differential equations arising from mean field type control can be 
performed with rather similar arguments as for MFG, see H for a work devoted to classical solutions. 
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The class of MFG with congestion effects was introduced and studied in [T^ in 2011, see also [ll[2] for 
some numerical simulations, to model situations in which the cost of displacement of the agents increases 
in the regions where the density is large. A striking fact is that in general, MFG with congestion cannot 
be cast into an optimal control problem driven by a partial differential equation, in contrast with simpler 
cases. In the present paper, we aim at studying mean field type control with congestion, in a setting in 
which classical solutions of the system of partial differential equations seem difficult to obtain. But, in 
contrast with MFG, mean field type control can genuinely be seen as a problem of optimal control of a 
partial differential equation. This will allow us to use techniques from the calculus of variations. Inspired 
by the works of Gardaliaguet et al, see we will introduce a pair of primal and dual optimization 

problems, leading to a suitable weak formulation of the system of partial differential equations for which 
there exists a unique solution. Note that is devoted to some optimal transportation problems (i.e. 
finding the geodesics for a class of distances between probability measures), whereas [6] deals with some 
special cases of MFG with possibly degenerate diffusions to which the above mentioned techniques from 
the calculus of variations can be applied. 


1.1 Model and assumptions 

This paper is devoted to the analysis of the second order system 
Ou 3H 

— (t, x) + vAu{t, x) + H{x, m{t, x), Du{t, x)) + m{t, x)— —(a:, m{t, x), Du{t, x)) = 0, (1.1) 

^(t,a:) - i^Am{t, x) + div(^m{t, ■), Du(t, ■))^ {x) = 0, (1.2) 

with the initial and terminal conditions 

m(0, a:) = mo(a:) and u{T,x) = ut(x). (1-3) 


Assumptions We now list the assumptions on the Hamiltonian il, the initial and terminal conditions 
Too £md ut- These conditions are supposed to hold in all what follows. 

HI The Hamiltonian H : T'^ x (0, +oo) x K is of the form 


H{x, rn,p) 


\P\^ 


■ £{x, m), 


(1.4) 


with 1 < /3 < 2 and 0 < a < 1, and where i is continuous cost function that will be discussed 
below. It is clear that H is concave with respect to p. Galling /?* the conjugate exponent of /?, i.e. 
/?* = /?/(/3 — 1), it is useful to note that 


H{x,m,p)= inf (^•p + L(a:,m,^)), (1.5) 

L{x,m,^) = {P — 1)P~^ +jt[x,m), (1.6) 

that L is convex with respect to and that 

L{x,m,^)= sup ■ p + H{x,m,p)). (1.7) 

peR"^ 


H2 (conditions on the cost i) The function £ : T'^ x K+ —K is continuous with respect to both variables 
and continuously differentiable with respect to m if m > 0. We also assume that m i—>■ m£{x^ m) is 
strictly convex, and that there exist q > 1 and two positive constants Ci and C 2 such that 


1 


1 


d£ , 

dm 



(1.8) 

+ C 2 . 

(1.9) 
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Moreover, since we can always add a constant to L, we can assume that 


i{x,m) >0, Va; £ Vm > 0. (1-10) 

The convexity assumption on m i—>■ m£{x,m) implies that m i—>■ mH{x,m,p) is strictly convex with 
respect to m. 

Moreover, we assume that there exists a constant C 3 > 0 such that 

\t{x,m) - (.{y,m)\ < + mf‘~^)\x - y\. (l-H) 


H3 We assume that P > q*. 

H4 (initial and terminal conditions) We assume that m-o is of class on T'^, that ut is of class on 
and that m-o > 0 and mo(x)dx = 1. 

H5 ly is a positive number. 

Remark 1.1. Note that Assumption [H5] can he relaxed: all what follows can be generalized to degenerate 
diffusions, i.e. to the following system of PDEs: 


du did 

— (t, x) + Aij {x) x) + H{x, m{t, x), Du{t, x)) + m{t, x) (x, m{t, x), Du{t, x)) 


dm, ^ 92 

dx.dxj 


/ dH \ 

{Aij{-)m){t,x) + dW[m{t, '^{t, •), Du{t, •)) j (x) 


0 , ( 1 . 12 ) 
0, (1.13) 


where {Aij)i<cij<d{x) = ^S(x)I]^(x) and Ti is a Lipschitz continous map from to with D 

possibly smaller than d. The necessary modifications can easily be found in m- 


1.2 A heuristic justification of (ll.l|l - (ll.3|l 

Consider a probability space (fl. A, V) and a filtration generated by a d-dimensional standard Wiener 
process {Wt) and the stochastic process (-^t)tG[o,T] in adapted to E* which solves the stochastic 
differential equation 

dXt = dt + dWt VtG[0,T], (1.14) 

given the initial state Xq which is a random variable J^'^-measurable whose probability density is mo. In 
(11.141) . is the control, which we take to be 


it = v{t,Xt), 


(1.15) 


where v{t,-) is a continuous function on T'^. As explained in [1], page 13, if the feedback function v is 
smooth enough, then the probability distribution mt of Xt has a density with respect to the Lebesgue 
measure, my{t, ■) £ PnL^(T'^) for all t, and my is solution of the Fokker-Planck equation 


dmy 


{t, x) — vXmy{t, x) + div(^my{t, pvlt, •)^ (x) 


= 0 , 


(1.16) 


for t £ (0,T] and x £ with the initial condition 

m„( 0 ,x) = mo(x), x £ 


We define the objective function 

rT 


J{v) = E 


f L{Xt,my{t,Xt),it)dt + UT{XT) 
Jo 


/ L{x,mv{t,x),v{t,x))mv{t,x)dxdt + / UT{x)my{T,x)dx. 

'[0,T]xT'i Jf^ 


(1.17) 


(1.18) 
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The goal is to minimize J(v) subject to (11.1611 and (11.1711 . Following A. Bensoussan, J. Frehse and P. 
Yam in [4], it can be seen that if there exists a smooth feedback function v* achieving J{v*) = min J7(t;) 
and such that rriv* > 0 then 

v*{t, x) = argmin^ {l{x, (t, x),v) + Vu{t, x) ■ w(t, a:)^ 

and (m„., u) solve dLHi, (HU and (lESli. The issue with the latter argument is that we do not know how 
to guarantee a priori that m will not vanish in some region of (0, T) x Hereafter, we propose a theory 
of weak solutions of (ITT]1-(IT311. in order to cope with the cases when m may vanish. 

Remark 1.2. Note that the system of partial differential equations that arises in a mean field game is 

Ou 

— ft, x) + vNult, x) + H(x, mit, x), Du(t, x)) = 0, 
at 

with im and To the best of our knowledge, for such a system with the Hamiltonian given in 

the existence of a solution is an open problem except in the stationary case with (3 = 2 and a = 1, see 
m; in the latter case, a very special trick can be used. Besides, the theory of weak solutions proposed 
below does not apply to MFG, because as explained above, MFG with congestion cannot be seen as an 
optimal control problem driven by a partial differential equation. 


2 Two optimization problems 

The first optimization is described as follows: consider the set ICq: 

ICo = {fGCmO,T]xT‘^):f{T,-)=UT} 

and the functional A on /Cq: 

= inf A{(j), m) 

m e T^((0,T) X T'^) 
m > 0 


( 2 . 1 ) 


where ^ 

A{4>,m) = J J m{t,x) ^-^{t,x) + i'Acj){t,x) + H{x,m{t,x), D(j){t,x))^ dxdt (2-2) 

+ / mo{x)(j){0, x)dx 
dT'i 

with the convention that if m = 0 then mH{x,m,p) = 0. Then the first problem consists of maximizing 

sup A{4>). (2.3) 

For the second optimization problem, we consider the set /Ci: 


(m,z) € i^((0,T) X T'^) X Li((0,T) x T'^; 

m > 0 a.e 

)Ci = { dm 

— - vAm + div2 = 0, 

at 

m(0, ■) = Too 

where the boundary value problem is satisfied in the sense of distributions. We also define 

mL{x,m,f^) if m > 0 

L{x,m,z) = ■^ 0 if {m,z) = (0,0) . 

+00 otherwise 


(2.4) 


(2.5) 
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Note that (m, z) !->■ L(x, m, z) is LSC on R x Using (II.6L (H2) and the results of [2] paragraph 3.2, 
it can be proved that (m, z) !->■ L{x, m, z) is convex on R x R^, because 0 < a < 1. It can also be checked 
that 

_ f svcp [—z ■ p + •mH{x,m,p)) if m > 0 or (m, z) = (0, 0) 

L{x,m,z)=l peR*^ (2.6) 

I +00 otherwise. 


Since L is bounded from below, / 

Jo Jfi 

(to, z) G /Cl. We are interested in minimizing 


L{x,m{t,x), z{t,x))dxdt is well defined in R U {+oo} for all 

(2.7) 


inf B{m,z), 


where if / / L{x,m{t,x), z{t,x))dxdt < +oo, 


JO 


B{m^z)= / / L{x,m(t,x), z{t,x))dxdt + / m{T,x)uT{x)dx, 


JO 


/jd 


and if not. 


B{m, z) = +00. 


( 2 . 8 ) 

(2.9) 


To give a meaning to the second integral in (12.8L we define w(t, x) = if rn{t, x) > 0 and w(t, x) = 

0 otherwise. From (HU) and (O, we see that L{x,m{t,x), z(t,x))dxdt < +oo implies that 

\w\^ G L^((0,T) X T"^), which implies that m\w\ G L^((0,T) x T"^). In that case, the 

boundary value problem in (12.4|) can be rewritten as follows: 


dm 


uAm + div(TOw) = 0, 


m( 0 ,-) = Too, 


( 2 . 10 ) 


and we can use the following Lemma which can be found in [6]: 

Lemma 2.1. If (to, w) G K-i is such that L{x, m(t, x), z(t, x))dxdt < +oo, then the map 1 1 —>■ m{t) 

for t G (0,r) and 1 1 — ^ toq for t < 0 is Holder continuous a.e. for the weak * topology o/P(T‘^). 

Remark 2.1. Following the proof of lemma 3.1 in m, we see that the Holder exponent in Lemma \2.1\ is 
greater than or equal to min • 

This lemma implies that the measure m{t) is defined for all t, so the second integral in (12.81) has a 
meaning. 


Lemma 2.2. 

sup .4(((i) = min B{m,z). (2-11) 

<^e/Co (m,2)e/Ci 

Moreover the latter minimum is achieved by a unique {m*,z*) G ICi, and to* G L‘^{{0,T) x T”*). 

Proof. Let us reformulate the optimization problem (12.31) : take Eq = C^([0, T] x T^) and Ei = C°([0, T] x 
T'^) X C°([0,r] X T‘^;R‘^). We define the functional T on Uq: 


= XTiJj) 



mo{x)(p{0, x)dx 


where xt{4>) = 0 if 4>\t=T 
A: Eo^ El hy 


ut and xt{4>) = +oo otherwise. Let us also define the linear operator 




For (a, b) € El, let Q(a, h) be defined by 

Q{a,b) = — inf f ( m{t,x) + H{x,m{t,x),b{t,x))) dxdt. (2.12) 

m G L^{{0,T) xT‘^) Jo 
m > 0 


Note that the infimum with respect to m is in fact a minimum: indeed, from (HH) and (HU) , we see 

that m I—>■ / / m{t,x) {a{t,x) + H{x,m{t,x),b{t,x))) dxdt is convex, coercive and continuous in the 

Jo Jt’^ 

set {m € L®((0,T) x : m > O}. Hence, this map is lower semi-continuous for the weak convergence 
in {m G L'^((0,T) x > 0}. On the other hand, since a minimizing sequence for ()2.12l) is 

bounded in L'^{{0,T) x T'^), we can extract a subsequence which converges weakly in L'^{{0,T) x T'^) to 
a nonnegative function. The weak limit achieves the minimum in (12.121) . We now aim at characterizing 
the optimal m. 

Let us first characterize 

K{x,'y,p) = mmip-y + pH{x, p,p)), (2.13) 

/i>0 

which is nonpositive and concave with respect to ( 7 ,p); since p i—?► pH(x, p,p) is C^, strictly convex on 
R+ and tends to -boo as p ^ -boo, we see that for any x G ii p ^ Q and 7 G R, or if p = 0 and 
7 -b £(x, 0 ) < 0 , then there exists a unique p = ip{x,j,p) >0 such that 

'y + H{x,p,p) + pHmix,p,p) = 0. 

Note that ifp = 0 and 7 -b^(a:, 0) < 0, then p = ^(x, 7 ,0) >0 is characterized by 'y+£{x, p)+p£m{x, p) = 0. 
We extend -0 by 0 in the set {( 7 ,0) : 7 + £{x, 0) > 0}. Therefore, 

K{x,'y,p) = il)[x,^,p)-f + ip{x,-i,p)H{x,ip{x,-i,p),p), (2.14) 

with the convention that mH(x, m,p) = 0 if m = 0. 

We claim that the map ( 7 ,p) !->■ ^(x, 7 ,p) is continuous in R x R'^. Indeed, 

1. the continuity of ( 7 ,m,p) 1 —?► j + H{x,m,p) + pHm{x,m,p) and the fact that this map is strictly 
increasing w.r.t. m in (0,-boo) implies that ( 7 ,p) 1 —>• ' 0 (x, 7 ,p) is continuous in R x R‘^\{0}. 

2. Similarly, the continuity of 7 1 —>■ 'ip{x,"f,0) stems from the continuity of the map ( 7 , to) i—> 7 -b 
i{x, m) -b mlrn{x, TO,) and its strictly increasing character w.r.t. m in (0, -boo). 

3. Let us prove that if (7,p) tends to (7,0) with p 0 , then 


lim '!/;(x, 7 ,p) ='0(x, 7 ,0). (2.15) 

(7,p)^(7,0) 

(a) If 7 -b (.{x, 0) < 0, then ^(x, 7 ,0) >0 and we get (I2.15|) from the same argument as in point 1. 

(b) Consider the case 7 + £(x,0) > 0. Suppose that ( 7 ,p) tends to ( 7 , 0 ) with p yb 0, and set 
p = 'ij}{x,^,p) > 0. We see that 

(1 — a)p~°'\p\^ = 7 -b w— (to !->• m£(x, m))(p) > 7 -b (.{x, 0 ) —>■ 7 -b £ix, 0 ). 
dm 

This implies that /2“|p|“^ is bounded as ( 7 ,p) —>■ ( 7 , 0 ), hence (I2.15|) . 

(c) Finally, we consider the case when 7 = —i{x,0) < 0 and ( 7 ,p) — >■ ( 7 , 0 ) with p 7 b 0; let us 
assume that for a subsequence, p is bounded away from 0 : passing to the limit in the identity 

(^(m to£(x,to))(/ 1 ) - ,^(x, 0 )) - {1 - a)p~°‘\pf = 7 - 7 , 


we obtain that 


lim -—(m m£(x,m))(p) = £(x,0), 

(7.p)^(7.o) dm 


which can happen only if p —7 0 and we reach a contradiction. Hence, (12.151) holds. 
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We have proved the continuity of ip with respect to The continuity of ip with respect to x follows 

from similar arguments, using the regularity assumptions on Therefore, ip is continuous in the set 

X R X 

It is also useful to notice that 


K(x,'y,p)= inf + H{x, p.p)) 

= inf ( M 7 + M inf (^ ■ P + L{x, p, ^)) 

^>0 \ 5 

= inf {pj + z ■ p + L(x, p, z 


and, from the Fenchel-Moreau theorem, see e.g. see m. that 


L{x,p,z)= sup {-p-f-z-p + K{x,-f,p)). 

(7,p)GRxR‘^ 


(2.16) 


Remark 2.2. Note that for all 0 p G and 7 £ R, the map p 1 —>■ p{j + H{x,p,p)) is strictly 
decreasing in some interval [ 0 ,p] where p> 0 depends on p and 7 , and that its derivative tends to —00 as 
p ^ O'^. Hence, if p* = 0 achieves the minimum of p^ pl^y + H{x, p,p)), then p must be 0. Similarly, 
7 must be such that 7 + i{x, 0 ) > 0 . 

With H2.14L the optimality conditions for (12.121) yield that 

Q[a,h) = — f f K{x,a(t,x),b{t,x))dxdt. (2-17) 

Jo Jt^ 

From (12.141) and the continuity of ip, we see that G is continuous on Ei. We observe that 


sup A{(p) = - inf {E{(p)+G{A{(p))) . (2-18) 

4ieKo 

By choosing (po{t, •) = wt, we see that iF{<po) < + 00 , G{A{(po)) < +00 and that t/ o A is continuous at (po- 
We can thus apply Fenchel-Rockafellar duality theorem, see [T5] : 

- inf {E{(p) + G{A{(p))) = min {E* {A* {m, z)) + G* {-m,-z)) (2-19) 

4>eEo {m,z)eE* 

where is the topological dual of Ei i.e. the set of Radon measures {m,z) on (0,T) x with values 
in R X R'^. If Eq is the dual space of Eq, the operator A* ; Fi* —^ E^ is the adjoint of A. The maps F* 
and G* are the Legendre-Fenchel conjugates of E and G- Following [5], we check that 


E*{A*{m,z)) 


J / UT{x)dm{T,x) 

( 9m . ,. 

if A — EAm + divz = 0, 

\ JtI 

[ m(0, •) = mo 

[ +00 

otherwise. 


where the boundary value problem is understood in the sense of distributions. 

On the other hand, from Rockafellar, [T5] Theorem 5, and (12.161) . see also [3, we see that 


G*{-m,-z) 



L{x, m^^{t, x),z^^{t, x))dxdt + 



/ 

V’ de 




where {m^^,z'^) and respectively denote the absolutely continuous and singular parts of 

(m, z), 6 is any measure with respect to which (to®'"®, is absolutely continuous, (for instance + 
and Loo{x, ■) is the recession function of L{x, •), i.e. 


Loo{x,m,z) 


sup —L{x, Am, \z) 
A>o A 


0 if {m,z) = ( 0 , 0 ), 

+00 otherwise. 
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Therefore, 


g*{-m, -z) = 


Hence, 


Jo 
+c» 


L(x, m{t, x), zit, x))dxdt 


if (m,z) e i^((0,T) X T'^) x L^{{Q,T) x 
otherwise. 

( 2 . 20 ) 


min {J^*{A*{m,z))+ g*{—m,—z)) = min B{m,z). 

{m,z)£E^ {m,z)GK.i 


and we obtain the desired result from (12.181) and (12.191) . 

Using (HSl), the strict convexity of m i—>■ m£{x,m) assumed in (H2), the convexity of the map (m > 
0, z) (see [2] paragraph 3.2), and the convexity of /Ci, we obtain the uniqueness of m* 

such that (m*, z*) G /Ci achieves a minimum of B for some z*. Moreover, from the strict convexity of 
{m,z) !->■ for to > 0, we deduce that ^ is unique in {{t,x) : m*{t,x) > 0}. Since z* = 0 in 

{(t^x) : m*{t,x) = 0}, the uniqueness of z* follows. It is clear from (H2) that to* € L‘^((0,T) x T'^). □ 


3 A priori estimates for a maximizing sequence of (l2.1l) 

Let M S M be the optimal value in m- Take a maximizing sequence {4>n)n^N for (EH). For some e > 0, 
it can be chosen in such a way that 

M -e< A{^n) < M. 

From the definition of M, we see that A{(f)n) < -4((/)„,too). Hence, 


M-e < 


J J rno{x) + vA(l)n{t,x) + H{x,mo{x),D(j)n{t,x))^ dxdt 

+ / mo{x)(f>n{0, x)dx 


10 Jt’' 


{—vDmo{x) ■ D(j}n{t, x) + mo{x)H{x, too(x), Dcj)n{t, x))) dxdt + / mo{x)uTix)dx. 


Ifd 


which implies that ||I2(/>ra||_L/9((o,r)xT<i) is bounded uniformly w.r.t. n. 
Let TOn achieve A{(j)n) = A{4>n,'mn)- Recall that to„ is unique. 

The optimality conditions for to„ in EH are 


+ vA(t>nit, x) + H{x, TO„(t, x),Dcj)n{t, x)) + mn(t, x)Hm{,{,x, m„(t, x),D4>n{t, x)) > 0 a.e., (3.1) 


0 = 


J J mn{t,x) + iyA(l)n{t,x)^ dxdt 


(3.2) 


JO JTd 
From (13.2L we deduce that 


TO„(t, x) {H{x, mn{t, x),D(j}n(t, x)) + TO„(t, x)Hjn{x, mn(t, x),D4>n(t, x))) dxdt. 


■^{4'n) = - / ml^{t,x)Hm{{x,mn{t,x),D(l)n{t,x))dxdt + / mo{x)(j)n{0,x)dx. 

Jo Jjd Jjd 


From (HI) and (H2), 


Hence, 


m'^Hjn{x,m,p) > —m'^ 

O 2 


C2m + m^-°‘\pf. 


[ [ (-^ml{t,x)-C2mn{t,x)+ml°‘{t,x)\D(l)n{t,x)f^dxdt-[ 

Jo JTd VO2 / Jjd 


< -M + e. 


mo(x)t/)n(0, x)dx 


(3.3) 









On the other hand, (13.IL (HI) and (H2) imply that for some constant C 4 > 0, 

+ C 4 m^~^(t,x) + C 4 > 0, a.e. (3.4) 

Multiplying (13.411 by for A large enough, integrating in (r, T) x T'^, we obtain that 

llL“(0.r;L<i*(T<i)) < C' (1 + ll”^"llL<i(\o,T)xT‘i) (T'')) ’ (^■^) 

and that 

'^^)IIl2((o,t)xT‘') ^ C' II IIL?* (T<i)) ' 

Remark 3.1. Note that the latter estimate does not hold with a degenerate diffusion as in RemMrk \l.l\. 
but it will not be used hereafter. 

This implies that 

J^^mo{x)(j)+{ 0 ,x)dx < C'(l|TO«|li;(\o-r),,TO + l) • (^.6) 

Combining (13.31) and (13.6L we obtain that 

I I (—m‘}^{t,x) — C2mn{t,x)+mlff°‘{t,x)\D(l)n{t,x)\^]dxdt+ [ mo{x)(j)~{ 0 ,x)dx 
Jo Jt’>' V^2 / Jjd 

< -M + e + C (l|mn|li<,(\o.T)xTO + 1) • 


The latter and (13.61) yield 

fl 

Jo Jt^ 

Let m(t,x) = lT<t- since A{(l)n,'rh) > A{ 4 >n), we obtain that 


mJfit.x) + m\ °‘{t,x)\D(l)n{t,x)\^) dxdt+ / \(j)n(f), x)\dx < C. 

JT<i 


(3.7) 


M -e < 


J J +i^A(^„(T,a;) + i7(a;, l,D(()„(T,a;))^ dxdt 

+ / mo{x) 4 >ni 0 ,x)dx 

JTd 

= / / H{xA, D4’n{T,x))dxdT + / (j)n{t,x)dx+ / (mo(a;) — 1)()>„(0, a;)c?x. 

Jq Jfi Jt** 

This implies that t 1 —>• 4>nit,x)dx is bounded from below uniformly w.r.t. n. Combining with (Ha, 

we get that 4 >~ is bounded in T°°( 0 , T; T^(T^)), and finally that </)„ is bounded in L°°( 0 , T; T^(T^)). 
Finally, setting = -£(-,to„) - TO„^(-,m„) and 0 < 7 ^ = (1 - l{m„>o}, we see that the 

sequence {'fn)n is bounded in L'? (( 0 ,r) x T"^), and that 


d4>n 

dt 


+ vA(j)n > 7n + In- 


Integrating the latter on (0,T) x T'*, we obtain that 


/ / ')^{t,x)dxdt < / UT{x)dx — / cj)n{0,x)dx— / / ^l^{t,x)dxdt, 

Jo Jt^ Jt^ Jt^- Jo Jt^- 


which implies that the sequence of positive function ( 7 ^)n is bounded in L^((0,T) x T'^). 
To summarize, we have proven the following lemma: 
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Lemma 3.1. The maximizinq sequence {(j)r,)nGn introduced at the heqinninq of^ \S.l\ is uniformly bounded 

Noting mn the nonnegative function achieving = A{4>n)y the sequence (m„)neN is uniformly 

bounded in L'^((0,T) x T'^). 

Calling ^n(t,x) = —H(x,mn(t,x),D(j)n(t,x)) — mnit,x)Hm{x,mn(t,x), D(j)n(t,x)), with the convention 
that H{x, m,p) + mHm{x, m,p) = £{x, 0) if m = 0 and p = 0, we can split 7„ as follows: 7™ = 7^ + 7^, 
where^l = to„) m„) andj^ = l{m„>o}- The sequence (7n)neN is uniformly 

bounded in L‘^*{{0,T) x T'^). The function 7^ is nonnegative and the sequence (7^)„gN is uniformly 
bounded in L^((0,T) x T"^). 


4 A relaxed problem 


4.1 Definition and first properties 

Let /C be the set of pairs (^, 7) such that 

• (t>&LP{Q,T-W^’^{T’^))CL°°{Q,T-M{T‘^)) 

• 7 G M{[0,T] X T'*) and 7^'^ = 7^ +7^, with 7^ G L'^ ((0,7) x T'*) is non positive, 7^ is non negative 
and 72 G Li((0,r) X , 7"“s G 7W+([0,r] x 




(4.1) 


It is clear that K, is convex. The following lemma implies that (j) has a trace in a very weak sense: 

Lemma 4.1. Consider (((',7) G 1C. For any Lipschitz continuous map ^ —>■ R, the map t 1— >■ 

Jjd ^ix)(p{t,x)dx has a BV representative on [0,T]. Moreover, if we note ^(x)<f)(t~^, x)dx its right 
limit at t € [0,7), then the map f, 1—>■ ^ix)(j){t'^,x)dx can be extended to a bounded linear form on 

C0(T‘^). 

Proof. Consider first a nonnegative and Lipschitz continuous function ^ ^ R+; the following identity 

holds in the sense of distributions: 

— f ^{x)(j){t,x)dxj — iz ( V(j){t,x) ■V£,{x)dx> f £,{x)^^‘^{t,x)dx, 

ht \ Jfd ) Jjd Jjd 


The second and third integral in the latter inequality belong to 7^((0,7)). From this, we deduce that 
1 1—>■ fj.a ^(x)(f(t, x)dx has a BV representative. 

If now ^ is a Lipschitz continuous function that may change sign, ^ —>■ M, then we write ^ 

and use the above argument separately for and we still obtain that t 1—>■ ^(x)if>(t, x)dx = 

fjj ^~^(x)(/>(t, x)dx — (x)cf>(t, x)dx has a BV representative. 

The continuity of ^ ^{x)(j){t'^, x)dx comes from the fact that (f G 7°°(0,7; AI(T‘*)). □ 

Thanks to Lemma |4.11 we may define the concave functional J on K, by 






K{x, 7^'^(<, x), 7)())(f, x))dxdt 



mQ{x)(j){0'^, x)dx, 


(4.2) 


and the relaxed optimization problem: 

sup J((^,7). (4.3) 

(0,7)gk: 

In (14.2L note that {x,t) >->• K{x,"f^^{t,x), D(j){t,x)) is a measurable nonpositive function, so the first 
integral is meaningful and has a value in [—oo,0]. 

Note also that, from (I2.14L it is possible to restrict ourselves to the pairs {(f, 7) G /C such that 7®''^(t, x) < 
—£(x,0) for almost every (t,x) such that D(j){t,x) = 0. Noting K, the set 


/C = {(0,7) G /C : 7’^'^(t, x) < —£{x, 0) for almost every (f, x) s.t. D(f){t, x) = 0} , 
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we have 


sup J{cj},j)= sup _ 7 ). 

((/),7)e/c (0,7)eK: 


(4.4) 


Lemma 4.2. For any (0, 7 ) G K, such that 


f j K{x,^^'^{t,x),D(j){t,x))dxdt > — 00 , 

Jo Jt’^ 


for any (to, z) G ICi such that m G L'^((0, T) x and 

i-T r 


to 


L{x,m(t,x), z{t,x))dxdt < + 00 , 


the following holds: 
for almost any t G (0,T), 


— / m{T,x)uT{x)dx + / m{t,x)<f>{t,x)dx + / / m{s,x) {^^^{sjx) + H{x,m{s,x),Dcj){s,x))) dxds 

JT<i JT-i Js=t 


< / L(x,m(s,x), z{s,x))dxds, 

J s=t Jt^ 


(4.5) 


/ m{t,x)(j){t,x)dx + / TOo(a;)((i(0''',a:)dx + / / in{s,x) {^^^{Sjx) + H{x,m{s,x),Dcj){s,x))) dxds 


< 


/ s =0 J 


L(x, to(s, x), 2 ( 5 , x))dxds, 


(4.6) 


and the meaning of fjj m(t, x)tf(t, x)dx will be explained in the proof. 
Moreover, if 


— / m{T,x)uT{x)dx + / TOo(x)(/)(0“'', x)(ia; + / / m{s,x) {"f^’^{s,x) + F[{x,m{s,x), D(j){s,x))) dxds 

Jt^ Jt>- j s=o Jt^ 

- / / L{x,m{s,x),z{s,x))dxds, 

J s=0 


(4.7) 


i/ien z{t,x) = m{t,x)Hp{x,m{t,x), D(l){t,x)) holds almost everywhere. 

Remark 4.1. Before proving Lemma \4.S\ note that the integrals of 

(s, x) I—to(s, x ) (7®'‘^(s, x ) + H(x, m{s, x), D(p(s, x))) in | 4 . 5 [ )- (|^. 7 | ) have a meaning in (—00, +00], because 
m(s,x) + H{x,m(s,x), D(j)(s,x))) > K(x,j^'^(s,x),D(j)(s,x)) a.e. and from the hypothesis of 

Lemma \4.‘d[ the latter is integrahle in (0,T) x with a nonpositive integral. 

Proof. We first extend to to [—1, T+ 1] by setting m{f) = toq for t < 0 and m{t) = m{T) for t > T. Note 
that m{T) is well defined from Lemma [2.11 Similarly, we extend z to [—1,T + 1] by setting z{t) = 0 if 

(o,r). 

Consider a regularizing kernel r]e{t, x) = e~'^~'^^i(^)/i 2 (f) where hi is a smooth even and nonnegative 
function supported in [— 1 / 2 , 1 / 2 ] such that hi(t)dt = 1, /12 is a smooth symmetric nonnegative function 
supported in [—1/2,1/2]'^ such that h 2 {x)dx = 1, and <5 will be chosen later. We define m 

and Ze = ? 7 e*z in (— 1/2,T + 1/2) x T'^. We can see that in (0,r) x T'^, 

- jzA (me Ite(o.T)) + divze = divR^ (4.8) 
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where = rj^-k [Dm lt^(o,T)) — ^tG(o,T)D{r]^-km). From [9], we know that —>■ 0 in F'*((0,T) x as 

e —>■ 0, because m G F'?((0,T) x T'^). 

We know that (l)(T) < ut- Thus (^{T, x)me{T, x)dx < UT(x)m^(T, x)dx and the latter integral 
converges to UT(x)m(T, x)dx in view of Lemma \TJ\ 

We know that mg —?► m in L'^{{0,T) x and that (j) G L^{{0,T) x T'^). From Assumption (H3), this 
implies that (ffm^ —>■ (f>m in L^((0,T) x T'*). 

From Lemma f4.ll we know that out of a countable set, 4’it~) = = 4’it) that x)m^[t, x)dx 

is well defined. 

From the latter two observations, up to the extraction of a subsequence, we may assume that 

/ cj){t, x)me{t, x)dx ^ / (j){t,x)m{t^x)dx 

Jfd. Jfd 

as e —>■ 0, with 4)[t~) = for almost all t G (0,r). Let t G (0,r) be such that the latter is 

true. 

Since G /C, 

T)'me{s, x))^ dxds 

— / UT{x)me{T,x)dx + / x)m^{t, x)dx (4-9) 

Jfi JT<i 

< — / / 'y^‘^{s,x)m^{s,x)dxds 

Jt JT^ 

On the other hand, from (14.8F 

I Id dxds 


It Jjd 

We deduce from (14.911 and (14.1011 that 


/ / D(f>{s,x) ■ Re{s,x)dxds + / / D(f>{s,x) ■ Ze{s,x)dxds 

Jt Jjd Jt Jfd 


(4.10) 


/ / D(j){s,x) ■ Re{s,x)dxds + / / D(f){s,x) ■ Ze{s,x)dxds 

Jt Jjd Jt JTd 

m^{s,x)H{x,m^{s,x),D(j){s,x))dsds — / UT{x)mi{T,x)dx + / (f>{t,x)me{t,x)dx (4.11) 

Jjd Jfd 


t Jfd 
T 


< - 


'y^'^{s,x)mi{s,x)dxds — / / m^{s,x)H{x,me{s,x),D(j){s,x))dsds. 

Jt Jjd 


h JTd 

Finally, we use the fact that 

i-T r 


h JTd 

pT 


{D(j){s, x) ■ Zc[s, x) — me{s, x)H{x, me{s, x), D(p{s, x))) dxds 


(4.12) 


> - 


Jt Jjd 

Combining this with (14.1111 . we find that 


L{x, me{s, x), Ze(s, x))dxds. 


— / UT{x)mi{T,x)dx + / (f){t,x)me{t,x)dx 

Jjd Jfd 


< - 


It Jjd 

pT 


me{s, x) (7®'‘^(s, x) + H{x, me{s, x), D(l){s, x))) dsds 


(4.13) 


+ / L{x,me{s,x), Ze{s,x))dxds + / / D(j){s,x) ■ Re{s,x)dxds. 


It JTd 


It Jjd 
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1. We have seen that the first line of (14.1311 tends to — UT(x)m(T, x)dx + x)m{t, x)dx. 


2. From Remark HU or more precisely the facts that K{x,j^'^{t,x), D(l){t,x))dxdt > —oo and 

that 

7 ®'‘^(s, x) + H{x, me{s, x),D(j){s, x))me{s, x) > K{x, x),D(j){s, x)), 
we can use Fatou lemma and get that 

/ / {"t^^{s,x) + H{x,m{s,x), D(j){s,x))) m{s,x)dxds 

Jt 


< lim inf 
e^O 


t JT'i 


( 7 ^'^(s,a;) + H{x,m^{s,x),D<j){s,x))) m^{s,x)dxds. 


3. D(l)(s,x) ■ Rt{s,x)dxds tends to 0 because tends to 0 in L'J((0,r) x T"^) and Dcj) € 
L'5((0,T) X T*^)) with P > q*. 

4. Finally, from the convexity of L with respect to (m, z) and from (ll.llF we see that 


lim sup 


e->0 Jt Jt>- 


L{x,me{s,x), Ze{s,x))dxds < 


It 


L{x, m{s,x), z(s, x))dxds. 


From (14.131) and all the points above, we deduce (14.51) . 
Similarly as for (14.1311 . we obtain that 


— / x)in^{t, x)dx + / (f>{0,x)me{0,x)dx 

Jjii JTi 

— / 'y^^{s,x)m^{s,x)dxds — / / m^{s,x)H{x,me{s,x),D(l){s,x))dxds 

Jo Jt<‘ Jo Jt^- 

ft 


(4.14) 


+ / / L{x,m^{s,x), z^{s,x))dxds + / / D(j){s,x) ■ Re{s,x)dxds, 


0 Jt^ 


0 Jt^ 


in which 6(0. x)mJn. x)dx has a meaninsf from Lemma 14.11 We claim that (^(0, x)me(0, x)dx 
Jja (j){0,x)mo{x)dx as e —?► 0: indeed, let ( £ (0,1) be the Holder exponent in Lemma o in view of 
Remarkwe know that C > min(l/2, (1 — a)//?); calling /i 2 ,e(x) = / 12 (f) and /ii,e(/) = e~^hi(-^), 

we get that for all x 


|m<;(0,a:) - (/i 2 ,e *mo)(a;)| = 


J jv) - wo(y))/ii.£(-s)/i 2 ,e(a; - y)dsdy 


<C J s^«-‘'-i/ii,£(-s)ds < 

Choosing S large enough, (i.e. such that — d — 1 > 0) and using the fact that toq is C^, 

lim ||/i 2 ,£*mo - mo||i=o(T<i) = 0. 
e^O 


Therefore, lime_>o ||me(0,-) — mo||ioo(Trf) = 0. The claim follows from the continuity stated in Lemma 

01 

The fact that x)m^(0, x)dx —>• fjj, </>(0, x)mo(x)dx and the arguments above imply (14.6F 

Let us now suppose that (14.71) holds: then the inequalities in (14.51) and (14.611 are equalities, for almost all 

t. For tr > 0 let us introduce the set 


E.(i) = 


(s,y)e[/,r]xT^ : 

m(s, y)ff(x, m(s, y),D(j){s, y)) - L{y, m{s, y),z{s, y)) < z{s, y) ■ D(l){s, y) - a 


If \E„{t)\ > 0, then for e > 0 small enough, > |^/(t(/)|/ 2, where 


(s,y)G[/,T]xT'' : 

'me{s,y)H{x,m^{s,y),D(j){s,y)) - L{y,meis,y), z^{s,y)) < Ze{s,y) ■ D(j){s,y) - f 
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Then (14.121) becomes 


It 

rT 


{Dif)[s, x) ■ Ze{s, x) — me(s, x)H(x, me{s, x), Dtp(s, x))) dxds 


(4.15) 


> - 


/ / L{x,me{s,x),Ze{s,x))dxds + -\Ecrit)\, 

It JT'i 4 


which implies that 


— / m{T,x)uT{x)dx + / m{t,x)<j){t,x)dx + / / m{s,x) {'-f’^'^{s,x) + H{x,m{s,x),Dcj){s,x))) dxds 

Jt^ Jf<>- Js=t 


< / L{x,m{s,x), z{s,x))dxds — —\E^{t)\, 

Js=t 4 

in contradiction with the fact that there is an equality in (ESI). 

Hence, m(s, y)H(x, m(s, y),D(j){s, y)) — L{y, m{s, y), z{s, y)) = z{s^ y) ■ D(j){s, y) holds almost everywhere. 
In view of 1)2.6L this shows that z{s,y) = m{s,y)Hp(x,m{s,y), D(j){s,y)) a.e. in {m > 0}. Furthermore, 


/ / L{x,m{t,x), z{t,x))dxdt < +oo, we see that z{s,y) = 0 a.e. in 

Jo Jt’>' 

{m = 0}. Hence, z{s, y) = m{s, y)E[p{x, m{s, y), D(j){s, y)) a.e. in (0, T) x T"^. □ 


from (1^ and the fact that 
{m = 0}. Hence, zi 

Proposition 4.1. 


sup = sup A{(l)). 

(0,7)G/C (pGK-o 


Proof. It is clear that \/(j) € Kq, —oo < A{(j)) < sup^^^.^jg^s^ J(^, 7 ): indeed consider m achieving A{(j), m) = 

A{(()), i.e. m = J) + vAtj), D(j^ . We know that < +00 and that m'^ < 00 . 

Let us take 7 = 7 ^ + 7 ^ where 7 ^ = —£(-,m) — and 7 ^ = ■^^^|ll<(>|^lm>o- It is easy to check 

that ((/), 7 ) € /C and that A{(j)) = J((/>,7)- Hence sup(^ J((/), 7 ) > sup^g^g^, A{(j)). 

For the reverse inequality, consider (</>, 7 ) G K such that J{<j),^) > — 00 . This implies that 

/ / K{x,^'^{s,x), D(j){s,x))dxds > — 00 . 

Jo Jt^ 

Let (m*, z*) be the pair of functions achieving ( 1 ^ . see Lemma O From Lemma El 

^ [ [ m*{t,x) {'y^'^{t,x) + E[{x,m*{t,x), D(l){t,x))) dxdt + ( mo{x)(j){Q'^,x)dx 

Jo dT-i 


< / , 

Js=0 

= B{m*,z*). 


L{x,m*{s,x),z*{s,x))dxds + / m*{T,x)uT{x)dx 


Hence, 


/jd 


min B{m,z)= sup (</>). 
(m,z)e/Ci (p^Ko 


and we conclude using (14.41) . □ 


4.2 Existence of a solution of the relaxed problem 

Proposition 4.2. The relaxed problem has at least a solution € /C. 

Proof. 
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Step 1 Consider the maximizing sequence 4'n for problem (1231) described in Lemma 13.11 and call 
the function such that A{(j)n) = A{4>n,'mn)- Also, let the functions 7 ^ and 7 ^ be defined as in Lemma 
13.11 We know that 

^ + vA(j)n =ln>ll,+ li, 

and {(j)n,Jn) S K. The definition of ensures that 7 ^ < —£(-,0). 

Up to the extraction of a subsequence, we may assume that ^ in L^(0, T, 1U^’^(T‘^)), ^ m in 

L«((0,r) X T'^), 71 ^ 7 ! inL«*((0,T) x T'^), with 7 I < -£(-,0) a.e., 7 ^ ^ 7 ^ in A^((0,T) x T'^) *. It is 
clear that 7 ^ is a positive Radon measure that we write 7 ^ = 7 ^ + 7 ®“s_ We see that 

d(j) . , 12 

— + vA(j) >7 = 7 + 7 . 

Thus for any function ^ in Lip(T‘^), the function t 1 —>• ^{x)(j){t,x)dx has a BV representative. More¬ 

over, since t !->■ Cn{t) = i{x)4>n{t,x)dx is bounded in L°°((0,T)) by C'||^||j;,oo(t<i), where C is in¬ 
dependent of we can assume (up to the extraction of a subsequence) that ^ C in A°°(( 0 ,T)) 
weak *, and |1 CIU“((o,t)) < C||?IIl°“(T‘^)- Since for all smooth function 1 /; : [ 0 ,T] —>• M, we know that 
^jd^{x)4in(iix)ip{t)dxdti&aAs to J,^d ^{x)(f>(t,x)'tjj{t)dxdt, we see that (^ = ^{x)4>{-,x)dx. There¬ 
fore \Jjd x)dx\ < C||^||ioo(Td). This shows that (^, 7 ) G 1C. 

Step 2 It can be proved that the map (^, 7 ) G /C K{x,—"f^^(t,x),—D4>(t,x))dxdt is the 

restriction to K, of the convex conjugate of the map 


G El 


lo JT<i 


L[x, x),z{t, x))dxdt. 


Therefore, the map (<(', 7 ) G /C !->■ K{x,^’^‘^{t, x), x))dxdt is upper semi-continuous for the 

weak * topology of L^{0,T; lUb^(T^)) x M, and therefore. 


K{x,j^‘^{t,x), D(l){t,x))dxdt > limsup 


Ar(x, 7 n(£, x), D(p{t, x))dxdt. 


(4.16) 


Step 3 Let = j.^d‘>xie{x)4>n{t,x)dx and C, = Jj.drno{x)(l){t,x)dx. We have seen that ^ C in 
L°°((0,T)) weak *; from the inequation satished by 


dC,n 

dt 


Dmo{x)D(l)n{t,x)dx - ( 7 ^(£, x) + 7 ^(£, a:))mo(x)dx > 0 , 

Jt’- 

which implies that 

-^(t) — V f Dmo{x)D(j)n(t,x)dx — ( 'y^{t,x)mo{x)dx > 0. 

dt Jrpd Jjd 

Thanks to the a priori bounds on and 7 ,), this implies that 

Ut)>Cn{o)-c{tw+t-^), V£G[o,r]. 

Letting n tend to + 00 , 

({t) > limsup Cn(0) — C{t~^ + £9 ), a.e. t G [0, T], 

n—¥oo 

Hence C(0) > limsup^^o^ Cn(0), i.e. J,j.d mo{x)4>{t, x)dx > limsup^^^ Jj.d mo{x)4>n{t, x)dx. 

Step 4 By combining the results of steps 2 and 3, we see that 

> limsup>l((/)n) = sup M(^) = sup J(^,7), 

ieKo (b.7)e/C 

Using the fact that 7 ^ < —£(-,0) a.e. and (j2.14ll . we can always decrease 7 ^ where Dcj) = 0 in such a way 
that ( 0 , 7 ) G 1C and the value of J{(j),^) is preserved. □ 
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5 Weak solutions to the system of PDEs 


Definition 5.1. A pair {4>,m) G L^(0,T;W^’^(T^)) x L‘?((0,T) x T'*) with m > 0 almost everywhere, is 
a weak solution of njfi-ms) */ 

1. m~°‘\D(f>\^lm>o € L^((0,T) X T*^), D(f>{t,x) = 0 a.e. in the region {it,x) : m{t,x) = 0}, 

e L^iiO,T) X T'^) and 

/ / m{t,x)L {x,m{t,x),Hp{x,m{t,x),D(j){t,x))) dxdt < + 00 . 

Jo Jt^ 

2. The following: 

-^{t, x) + i^A(j)(t, x) + H(x, mit, x), D(j)(t, x)) + mit, x)Hm(x, m(t, x), Dcfit, x)) > 0, 
ot 

(l>(T,x) < urix), 

holds in the sense of distributions, with the convention that 

H{x,m,p) + mHmix,m,p) = i{x,0) if m = 0 and p = 0. 

3. The following: 

dm, 
dt 

m(0, x) = mo{x) 
holds in the sense of distributions 

4- There exists a constant C such that for each function ^ G Lip{T‘^), 


-ft, x) — vAmft, x) + div(m(t, ■)Hpf,mft, ■),D(j)ft, •)))(a;) = 0, 


II f ^{x)(f{t,x)dx\\L^(^o,T) 

JT<i 


5. The following identity holds 

fT r 


(5.1) 


(5.2) 


(5.3) 


0 = 


/o 


mft,x)L(x,mft,x),Hp[x,mft,x),D(j)ft,x)))dxdt + / / m^ft,x)Hm{x,mft,x),D4>ft,x))dxdt 


/O 


+ / m{T,x)uT{x)dx — / mo{x)(l){Q'^ ,x)dx 


(5.4) 


Theorem 5.1. Let fm*,z*) G K-i be the minimizer of and (</)*, 7 *) G 1C be a maximizer of 

i4L^ . Then is a weak solution of il. D - HTB) . and z*ft,x) = m* ft, x)Hp{x,m* {t, x), Dcf)* ft, x)), 

'y*'^‘^ft,x) > —H{x,m*{t,x),D(f*ft,x)) — m*ft,x)Hm{x,m*ft,x),D(f>*ft,x)) almost everywhere. 
Conversely, any weak solution {(fjfh) of SI. 1\) - S1.3\) is such that the pair ffh,rhHp(x,rh, Dcf)) is the min¬ 
imizer of and (</), 7 ) is a maximizer of with 

fft, x) = —H{x, mft, x),D(j)ft, x)) — fhft, x)Hm{x, fhft, x),D(i)ft, x)). 


with the convention (EJP. 

Proof. Let fm*,z*) G /Ci be the minimizer of (12.711 . (which implies that m* G L‘>{{0,T) x T*^) and that 
Jq Jjd L{x,m*{t,x), z* ft, x))dxdt < + 00 ) and (<()*, 7 *) € /C be a maximizer of (14.31) . We know that 
Jo Ijd K{x,^*’^^ ft, x), Def* {t,x))dxdt > — 00 . Hence, 


/O dT'i 
rT 


L{x,m*ft,x),z*ft,x))dxdt-\- / m*fT,x)uT{x)dx 

Jfi 


f j K{x,^*'’^{t,x),D(j)* ft,x))dxdt -\- f mo{x)(l)*{0,x)dx 
Jo Jt^ JT‘* 


< / / m*{t,x) {'y*’^^{t,x) T[{x,m*{t,x),D(j)*{t,x))) dxdt / mo(x)(()*(0, a:)dx 

Jo Jt’>' dT-i 
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But from Lemma 14.21 and especially (14311-(1461). we see that the latter inequality is in fact an equality; 
then, from the last part of Lemma l42l this yields that z*{t,x) = m*{t,x)Hp{x,m*{t,x),D(j)*{t,x)) and 
that 


L{x, m* {t, x), z*{t, x)) 


L {x, m*{t, x),m*{t, x)Hp{x, m*{t, x), x)) 


almost everywhere. We have proved that m*Hp{-,m*, Dcj)*) S L^((0,T) x and that (15.311 holds. 
Moreover 

pT p rT r 


/ / K{x,^*’^{t,x),D(j)* {t,x))dxdt 


m*{t, x) x) + H{x, m*{t, x), x))) dxdt. 


Jo JT'i 


Jo Jv 


This implies that 


K{x,j*'^‘^{t, x),Dcj)*{t,x)) = m*{t, x) ( 7 *’^'^(t,x) -|- H{x, m*{t,x), D(j)*{t,x))) (5.5) 

almost everywhere. In view of Remark 12.21 this implies that Dcjfit^x) = 0 almost everywhere in the 
region where m*(t, x) = 0 , and that 

'y*’‘^^{t,x) =—H{x,m*{t,x),D(j)*{t,x)) — m*(t,x)Hmix,m*{t,x),D(j}*{t,x)), a.e. where >0 

and with the convention (15.2L we see that 

x) > —H{x, m*{t, x), D(j)*(t, x)) — m*(t, x)Hm{x, m*{t, x),D4)*(t, x)), a.e.. 

These observation imply that S L^((0,T) x T'^). Using also (15.511 . we see that 

€ L^((0,T) X T'^). From the inequalities satished by 7 *’'^'^, we see that (15.111 holds. 
Finally, dSH) is obtained by combining all the points above. We have proved that (</>*,m*) is a weak 
solution in the sense of Definition ICT 


Suppose now that (</>, m) is a weak solution in the sense of Definition 15.11 Let us choose 'y(t,x) = 
—H{x,fh(t,x),D(l){t,x)) — m(t,x)Hmix,m{t,x), D<l){t,x)) always with the convention (15.2p . (we split 7 

as follows: 7 = 7 ^ -p 7 ^, where 7 ^ = 7 ^ = (1 - a)^-^^l{fh>o}) and z(t,x) = 

fh{t,x)Hp(x^m{t,x), D(j)(t,x)). It is clear that (^,7) € 1 C and that (rh,z) € ICi. 


From the definition of 7, we see that a.e., m{t,x) = argmin^>Q (//7(t, x)-I-/iFl(x,/i, D(/)(t, x))) . 
Moreover, since a.e., Hp{x, fh{t, x),D(l){t, x)) = argmin^ (^ • x) + L(x, fh(t, x), ^)), we see that a.e. 

(m(t, x),z(t, x)) = argmin(^ x) + z ■ x) + Z(x, /r, z)^ , i.e. 

K{x, 7(t, x),D(l){t, x)) = fh(t, x)7(t, x) -I- z{t, x) ■ D(l){t, x) -I- L(x, m(t, x), z{t, x)). ( 5 - 6 ) 


Note also that 0 < fj' 1x1(1, x) 7 ^(t, x)dxdt < -l-oo and that 


f ( \z{t,x) ■ D(l){t,x)\dxdt < C f f °'{t,x)\D(j){t,x)\^dxdt < +00. 
Jo JTi JO JTi 


10 JTi JO JT'i 

Let (m* ,z*) be the solution of (13711 . From Lemma m 


B(m*,z*) = 


10 JTi 

pT 


L(x,m*(t,x),z*(t,x))dxdt + / m* (T,x)uT(x)dx 


Ifd 


> 


/ m*{t,x) (^(tjx) + H(x,m*(t,x),D(j){t,x))) dxdt + / mo{x)(j){Q,x)dx. 
Jfl JTi 


Hence, 


B{m*,z*)> I I K(x,^(t,x), D(j)(t,x))dxdt + J mo(x)(j)( 0 ,x)dx 


lo JTi 

lo JTi 


m(t,x) ('-f{t,x) + H{x,m(t,x),Dcj){t,x))) dxdt + / mo(x)(()(0, x)dx. 
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where the first line comes from (15.611 and the definition of K, and the last line comes from the definition 
of 7 . 

Finally, from the latter inequality and (15.4L we deduce that 


B{m*,z*)> f [ L{x,m{t,x),z{t,x))dxdt 

Jo 

=B{fh, z). 



m{T, x)uT{x)dx 


Therefore (m, z) achieves the minimum in (EJl) and m = m*, z = z*. 

It remains to prove that (<^, 7 ) achieves the maximum in (14.31) . We deduce from Lemma [521 Proposition 
O and the latter point that 


max J{(j),j)=B{m,z) = 
(0,7)e/c 


/ L{x,m{t,x),z{t,x))dxdt + / m{T,x)uT{x)dz 
Jfii 


/o 


m^{t, x)Hm{x, m{t, x), D(j){t, x))dxdt 


Ifd 


x)dx, 


where the last line comes from (15.41) and the definition of z. But, using the definition of 7 , 


—m^(f, x)Hjn{x, m{t, x), D(j){t, x)) =m{t, x)^{t, x) + m{t, x)H{x, m{t, x),D(j){t, x)) 

=m(t, x)^{t, x) + z(t, x) • D(j)(t, x) + L(x, fh{t, x), z{t, x)) 
=K{x,j{t,x),D^{t,x)) 


where the second line is obtained using the definition of z and the third line is obtained using (15.611 . 
Combining the latter two observations, we see that 

nT r r 


max J((/), 7 )= / / K{x,j{t,x),D(j){t,x))dxdt+ / mo{x)(f){0'^,x)dx = 

Jo Jt'^ Jt^ 


which concludes the proof. □ 


Theorem 5.2. There exists a unique weak solution of 

Proof. Existence is a direct consequence of Theorem 15.11 Proposition [42] and Lemma [221 

From Theorem 15.11 we also see that any weak solution (0 ,to) of (I1.1D - (I1.3D is such that the pair 

{m,mHp{x,rh, D(f)) is the minimizer of (12.711 . Thus {fh,fhHp{x,rh, Dcf)) is unique. 

Consider now two weak solutions {(j)i,fh) and ( 02 , th) of (ll.lll - (ll.3|l . We know that mHp{x,fh, Df)i) = 
rfiHp{x,rfi,D(j) 2 )'. this implies that |D 0 i(t,x)|^“^D 0 i(t,x) = \D(j) 2 {t,x)\^~'^D(j) 2 {t,x) at almost every 
{t,x) such that fh{t,x) > 0 : hence at almost every {t,x) such that m{t,x) > 0 , |D 0 i(t,x)| = |D 02 (t,x)| 
and finally D0i(f,x) = DJ> 2 {t,x). Moreover, D0i(f,x) = D(j) 2 {t,x) = 0 at almost every {t,x) such that 
fh{t, x) = 0. Therefore, Defift, x) = D 02 (i, x) at almost every (t, x). This means that 0i —02 only depends 
on t, and that 71 = 72 a.e., where 7 i(t, x) = —H{x, ffilt, x),D(j)i(t, a;)) — fh{t, x)Hm{x, fh(t, x), D(j)i{t, x)), 
using convention (EH). 

Using (El, we deduce from the previous points that 0i(t = 0"*", •) = 02 (t = O'*", •) a.e. in T^. 

We set 7 = 7 i = 72 S L^((0,T) x T'^) and z = fhHp{x,m, Df>i) = fhHp{x,rh, D4>2). Going back to the 
proof of Lemma 1321 we see that for almost all t, both ( 0 i, 7 ) and ( 02 , 7 ) achieve the equality in (14.61) 
with (m,z) = ( 771 , z). From the previous points, this implies that 



mit, x)4>i{t, x)dx 


Ijd 


mft, x)(j) 2 (t, x)dx 


for almost every t. Since 0i(t,a;) — 02 (t,a;) only depends on t and m{t) is a probability measure, the 
latter implies that 0 i = 02 holds almost everywhere. □ 
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